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We develop a method to entangle neutral atoms using 
cold controlled collisions. We analyze this method in two 
particular set-ups: optical lattices and magnetic micro-traps. 
Both offer the possibility of performing certain multi-particle 
operations in parallel. Using this fact, we show how to im- 
plement efflcient quantum error correction and schemes for 
fault-tolerant computing. 



I. INTRODUCTION 

Entanglement is one of the most intriguing features 
of Quantum Mechanics. However, there are very few 
physical systems in which entanglement can be system- 
atically studied in a controlled way. Those systems in- 
clude ion-traps cavity QED photons (H-pl, 
and molecules in the context of NMR [p6|-29[ (see |3C| ] 
however). Very recently, we have identified a new way 
of entangling particles by using cold controlled collisions 
with which one could study ex per imentally basic issues of 
quantum information theory [pl| . Given the impressive 
experimental advances made so far in the fields of neu- 
tral atom trapping and cooling 
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of Bose Einstein condensation (BEG) of ultracold gases 
p6|-p| , that proposal opens a new perspective to several 
experimental groups who so far have concentrated their 
efforts in other fields of Atomic Physics. 

In the present paper, we build upon the work in |3|] 
and explore the idea of using atomic controlled cold colli- 
sions for entangling neutral atoms in optical lattices (see 
also [^) and in arrays of magnetic micro-traps. We show 
how to perform two-qubit gate operations with those sys- 
tems obtaining very high fidelities. We propose a variety 
of experiments to entangle particles using state-of-the- 
art technology. We also concentrate on the unique possi- 
bilities that these set-ups offer to perform multi-particle 
entanglement operations in parallel p^-]4^]. Using such 
parallelism, we show how to implement efficient error cor- 
541 and fault-tolerant quantum computation 

% 

The paper is organized as follows. In Sec. II we dis- 
cuss the use of ultracold collisions as a mechanism for 
entangling neutral atoms. Such collisions can be brought 
about by either moving the potentials in certain spatial 
directions or by modifying the shape of the trapping po- 
tentials. In Sec. Ill we describe two systems in which 
such operations can be implemented. These are optical 
lattices p2|,|63|-p6l and magnetic microtraps ||67]-|7l| both 



rection 
schemes 




of which have been studied experimentally in detail in 
the past. In Sec. IV we describe a class of multi-particle 
entanglement operations that can be realized in these 
systems (we concentrate here on optical lattices). The 
usefulness of such operations for quantum computing de- 
pends on certain conditions that need to be satisfied in 
an experiment. Among these conditions, the filling prob- 
lem, i.e. how to fill the potentials with regular patterns 
of atoms, is most outstanding. We discuss these matters 
and show that even under present-day experimental con- 
ditions, very interesting entanglement studies could be 
performed. Section V summarizes the main results and 
discusses their relevance for future research. 



II. ENTANGLEMENT OF ATOMS VIA COLD 
CONTROLLED COLLISIONS 

In this Section, we consider two bosonic neutral atoms 
with two internal states trapped by conservative poten- 
tials and cooled to the motional ground states. Initially 
these two particles are sufficiently far apart so that they 
do not interact with each other. We then assume the 
shape of the potentials to be varied in a way that depends 
on the internal state of the atoms so that the two parti- 
cles come close to each other if they are in certain internal 
states. As we will show, this can be done e.g. by moving 
the center position of the trapping potentials state se- 
lectively, or by switching off a potential barrier between 
the two atoms for one of the two internal states. In both 
cases the particles will interact via s-wave scattering with 
each other in a coherent way when they are close to each 
other. After the interaction has taken place the particles 
are restored to their initial position. In this way one can 
implement conditional dynamics and realize a fundamen- 
tal two-qubit gate. 

Note that we are dealing with bosons. Therefore, we 
have to use symmetrized wave functions for describing 
the two particles. It will turn out that if the center po- 
sitions of the trapping potentials are moved state selec- 
tively, particles in the same internal state will always 
be so far apart that their wave functions never over- 
lap. Thus, we will not care about the symmetrization 
in this case. On the other hand, if the potential bar- 
rier is switched off for one internal state, particles in the 
same internal state will come close to each other and 
symmetrizing the wave function is essential. 
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A. Hamiltonian 

Here we deal with the interaction Hamiltonian of two 
neutral atoms 1 and 2 with internal states |a)i_2 and |6)i,2 
trapped by conservative potentials V^^ (x^, t) whose func- 
tional dependence on the coordinate x^, with i = 1,2 the 
particle index, depends on the internal state of the par- 
ticle Pi^2 = o,,b. Initially, the two particles are in the 
ground state of the trapping potentials and the centers 
of the two potential wells are sufficiently far apart so that 
the particles do not interact. Then the form of the po- 
tential wells is changed such that there is some overlap 
of the wave functions of the two atoms, and the particles 
will interact with each other. This interaction between 
the atoms in two given internal states Pi and (32 can be 
described by a contact potential 



(xi - X2) 



5'(X1-X2), 



(1) 



where afi'^^ is the s-wave scattering length for the corre- 
sponding internal states describing elastic collisions and 
m is the mass of the particles. This zero energy s-wave 
scattering approximation will be valid as long as we as- 
sume that Uosc , the rms velocity of the atoms in the vibra- 
tional ground state, approximately given by Uqsc ~ ^0^1 
is sufficiently small . Here an is the size of the ground 
state of the trap potential, and lu is the first excitation 
frequency. Thus we can describe the evolution of the 
system by the Hamiltonian 



H ^ J2 ^^^^^ ® \(^' 

01,02 
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(/?i|® 1/32)2 (/92 



(2) 



where 

fJ01p2 ^ ^ 



2m 



u''i*(xi - X2). 



(3) 



Here is the momentum operator. 



1. Interaction in perturbation theory 

We want to treat the interaction term in the Hamilto- 
nian Eq. (^) perturbatively. For particles in two different 
internal states /3i ^ /32 we find the correction to the en- 
ergy due to the interaction as 



m 



dxHl^* (x,i) 



(4) 



where V'f ' (x, t) is the normalized one-particle wave func- 
tion of particle i in internal state (3i in the time de- 
pendent potential V^^ (x, t) . If the particles are in the 
same internal state (3i = P2 — P, we have to account for 



the Bose statistics i.e. use the properly normalized sym- 
metrized two-particle wave function for calculating the 
energy shift. We therefore find 



m(l + lap) 



dx[] hAf (x,t) 



where 



— dx 



(X,t))*^2^(x,i), 



(5) 



(6) 



For general /3i , P2 we find the phase accumulated due to 
the interaction in the time interval [— r, r] by 



B. Moving potentials 



(7) 



One way of controlling the interaction between the par- 
ticles is to move the center position of the potentials 

towards each other in a 



state-dependent way while leaving the shape of the po- 
tential unchanged. By moving the potential we get two 
kinds of phase shifts. A kinetic phase which is a single- 
particle phase due to the kinetic energy of the particles 
and an interaction phase due to coherent interactions be- 
tween two atoms. First we will define these two phases 
for general trapping potentials and afterwards specialize 
them to moving harmonic potentials. Finally, we will 
show how conditional dynamics can be realized. 



1. Kinetic phase 

First we want to consider a single atom in internal 
state 1/3) trapped in the instantaneous ground state ipQ 
of a moving potential well y(x — x''(t)). The center po- 
sition of the potential is moved along a trajectory x^(i). 
Ideally, we want the atom to remain in the ground state 
of its trapping potential and to preserve its internal state 
during the motion. This corresponds to the transforma- 
tion from t — — r to t = r 



i^o[x-x^(-r)]->e-*^'V^o[x-x'3(r) 



(8) 



where the atom remains in the ground state of the trap- 
ping potential and preserves its internal state. Trans- 
formation (H) can be realized in the adiabatic limit [ [73[ , 
where we move the potentials so that the atoms remain 
in the instantaneous motional ground state. Adiabaticity 
requires |x^(i)| <C Wosc for all times t. The phase cj)^ can 
be easily calculated in the limit |x'^(i)l Vosc/t. We 
find the kinetic phase 



m 
2h 



dt 



(9) 
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2. Interaction phase 

Let us now consider two particles i — 1 , 2 in differ- 
ent internal states trapped in the ground states of 
two moving potentials. Initially, at time t = — r, these 
wells are centered at positions x^, sufficiently far apart 
(distance ci = Xi — X2) so that the particles do not in- 
teract. The positions of the potentials are moved along 
trajectories xf ' (t) so that the wave packets of the atoms 
overlap for certain time, until finally they are restored 
to the initial position at the time t = t. We assume 

that: (i) |xf'(t)| ^ Vosc (adiabatic condition) so that the 
particles remain in the ground states of the moving trap- 
ping potentials; (ii) The interaction can be treated per- 
turbatively, where \AE^^^^{t) \ <C hcj so that no sloshing 
motion is excited. In that case, we realize the transfor- 
mation 



'0O(X1 - Xi)V'o(x2 - X2) ^ 

e~"^'0o(xi - xi)V'o(x2 



where = cf)'^ + cf)'^ + 
0/31/32 defined in Eq. (0). 



X2), (10) 
with the collisional phase 



Xj -|- (5x^'(t). As shown in Fig. Qthe trajectories 5yi^^{t) 
are chosen in such a way that 5:k!^^{—t) — 5x^^{t) — 
and the first atom collides with the second one only 
if they are in states \a) and respectively (|xj(i) — 
X2(i)| ^ ao Vi). This choice is moti vated by the physical 

The fact that 



implementation considered in Sec. Ill A 
Xi does not depend on the internal atomic state and the 
shape of the two potentials is the same at times ±r allows 
one to easily change the internal state at times t = ±r by 
applying laser pulses. If the conditions stated above are 
fulfilled, depending on the initial internal atomic states 
we have: 

\a)i\a)2^e~^'*''\a),\a)2, 

|a)i|6)2-e-(*°+*'+*°'')|a)i|6)2, 

|&)i|a)2->e-*(^"+^')|6)i|a)2, 



\b)i\b), 



\b)i\bh 



(12) 



3. Moving harmonic potentials 



where the motional states remain unchanged. The kinetic 
phases (p^ and the collisional phase 4'°'^ can be calculated 
as stated above. We emphasize that the 0'^ are (trivial) 
one-particle phases that, as long as they are known, can 
always be incorporated in the definition of the states \a) 
and 1 6). This realizes a fundamental two-qubit quantum 
gate for certain values of 0"'', e.g. = tt. 



Here we specialize to harmonic trapping potentials. 
The wave function i/jf ' (x, t) of a particle in a moving 
harmonic potential can be found analytically. In the 
Appendix ^ we show that when we start to move the 
harmonic potential at time — r with the particle in its 
motional ground state and stop to move the potential at 
time T, the condition for the particle to end up in the 
motional ground state at r is given by 



(11) 



This condition is weaker than the condition |xf'(f)| ^ 
Wosc for adiabaticity, and means that the particle need 
not be in the instantaneous ground state of the moving 
potential at all times, but only at the final time. The 
kinetic phases can be found exactly (cf. Eq. (^). If 
\AE^^^^ {t)\ <C Uuj is satisfied, the interaction phase can 
be found by Eq. ^ since the ipf' (x, t) are known. It 
is also possible to generalize these results to the case in 
which the trap frequency changes with time . 



4- Implementation of conditional dynamics 

Let us now assume that we can design the potentials 
such that atoms in the internal state experience a 
potential V^'{xi,t) — V{xi — xf'(i)) which is initially 
-t) centered at position x^. We assume that we 



it 

can move the centers of the potentials as follows: xf ' {t) — 




xi(i) 



FIG. 1. Configurations at times ±r (a) and at t (b). The 
solid (dashed) curves show the potentials for particles in the 
internal state \a) {\b)). Center positions xf'(f) and displace- 
ments Jx''' {t) as defined in the text. 



C. Switching potentials 

The interaction between the particles can be controlled 
also in another way, for example by changing with time 
the shape of the potentials depending on the particles' in- 
ternal states. Different regimes for the time-dependence 
of the potential are possible. The two limits of ex- 
tremely slow (adiabatic) or extremely fast (sudden) po- 
tential changes are both interesting and lead to peculiar 
schemes. Here we will analyze the latter case. We con- 
sider two atoms initially trapped in two displaced wells. 
At a certain time the barrier between the wells is sud- 
denly removed in a selective way for atoms in state |6), 
whereas it remains unchanged for atoms in state \a). The 
atoms are allowed to oscillate for some time, and then the 
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barrier is raised again suddenly such as to trap them back 
at the original positions. During the process they will 
acquire both a kinematic phase due to the oscillations 
within their respective wells, and an interaction phase 
due to the collision. We will calculate such quantities 
and look for the optimal switching time required in or- 
der to maximize the fidelity for a quantum gate relying 
on this scheme, which we will estimate quantitatively for 
the relevant physical example in Sec. IIIB. 



1. Kinematic phase 

Let us first consider the time-independent problem of 
an atom subject to a three-dimensional potential whose 
functional form along x depends on the internal atomic 
state (3 ~ a,b: 



V'^{^)^v'^{x)+v^{y) + v^{z). 



(13) 



Here the u's are single-well trapping potentials, and 

are centered around x — xq, a; = 0, respectively. We 
assume that the atom is initially prepared in the motional 
state 1^'+), where 



1'+(x) = V+(a;)V'i(y)V'±(^) 



(14) 



and ip+j V'J- are the ground-state wave functions of w", 
?;_L with eigenvalues E"", E± respectively. Thus ^'-i-(x) 
is peaked around the position Xq = (a:o,0,0), coinciding 
with the center of y°(x) but displaced from the one of 
y''(x). Therefore, if the atom is in internal state |a), its 
motional state after a time t will be unchanged up to a 
phase 0° = (i?° -I- 2E±)t/h. If it is instead in state \b), 
it will start oscillating within the well, thus picking up a 
different phase (j)^ due to the kincmatical evolution, and 
possibly coming back at the initial position after some 
time. 



replacing the Hamiltonian (||) in Eq. Here w^{x,t) 
is a combination of the v^{x) whose form changes with 
time, and u^'-^^ is an effective interaction potential tak- 
ing into account the integration over y and z, and there- 
fore depending on the shape of v±. We shall study the 
dynamics at t > for different values of (/3i,/32) sepa- 
rately, li /3i = /32 = P the total initial normalized state, 
symmetric under particle interchange, is 



|V;+)i|V^_)2 + |V^-)i|V^+)s 

V2 



(16) 



where the initial overlap (-0- 1 has been neglected 
in computing the normalization. If both particles are in 
state I a), no interaction takes place and thus the colli- 
sional phase 0"° — 0. Therefore, we shall now consider 
in more detail the situation in which both particles are 
in state \b) and thus move within the well v^. In the 
absence of interaction, after an oscillation period T they 
would come back exactly to the initial state. Due to the 
interaction, two effects arise: an additional phase, which 
is accumulated by the wave function as the number of 
undergone oscillations increases; and a slight decrease in 
their frequency, because the atoms acquire a small delay 
in their motion inside the trap as they come out from 
a collision. These effects have to be evaluated in detail, 
since they influence the attainable fidelity for a quantum 
gate based on this scheme. For symmetry reasons, the 
relative coordinate motion decouples from the center of 
mass motion, which is not affected by the interaction and 
can be solved analytically. For an explicit calculation, it 
is now needed to specify the form of the potentials in 
Eq. P). 



3. Switching harmonic potentials 



2. Interaction phase 



In order to perform the calculations analytically, the 
potentials in Eq. (hsh are chosen to be harmonic: 



We now consider two atoms 1 and 2 initially (at t = 0) 
prepared in the motional states |^'+) and l^*-), the latter 
being defined as in Eq. ( |l^ but with ip- (x) = ip+ {—x) re- 
placing ijj+ix)- We assume that the particles are subject 
to the potentials J2<;=+ -'^i'^^W^'i'^^i): where Q de- 
notes the step function. If any one of them is in state \b), 
for t > it will start oscillating within the well, eventu- 
ally interacting with the other one. If v± is much steeper 
than , then the probability of transversal excitations 
can be neglected, i.e. each atom remains in the ground 
state along y and z. By integrating over these variables, 
the problem is then reduced to a one-dimensional two- 
particle Schrodinger equation, with 



n 



0102 



E 

i=l L 



2m 



W^'{Xi,t) 



+ U^^^^{xi-X2) (15) 



v^x)='-^{x-xor, 
v\x) = — -a;2. 



(17a) 
(17b) 
(17c) 



where lo± ^ > oj. Our scheme for gate operation is as 
follows: initially the two particles are separately stored in 
two displaced harmonic wells at ±a;o as described above, 
i.e. with the potential (Fig. ||a) 



u.''(a;,t < 0) = ^ Q{^x)v''{qx), 



?=+,- 



w^{x,t< 0) = w°'{x,t < 0) 



(18a) 
(18b) 
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in the one-dimensional Hamiltonian Eq. ([Tsl). At t = 
the potential undergoes a sudden change, namely the 
barrier between the two wells is selectively switched off 
for state \b) only (Fig. ^3): 



^"(2;, < t < t) = w^ix, t < 0); 



(19a) 
(19b) 



Then at t = t, the potential barrier is suddenly restored: 
w°-'^{x,t > r) = w°-'''{x,t < 0). The time evolution at 
< t < r is characterized by oscillations with periodic- 
ity T = 2tt/uj. The projection of the evolved CM wave 
function on the initial one 



K^^''MwiV'^'M(o))r= 



■ sin^ (tot) 



(20) 



has instead a period of T/2, because of the parity of the 
spatial wave function. The time-dependent energy shift 
(0) due to the interaction turns out to be 



AE''''{t) = afhu 



bbi 



(21) 



where f2(t) = lj^ ujq / [uj'^ cos^ {ujt) + LjQsm'^{ujt)]. Hence 
the interaction-induced phase shift accumulated after 
each oscillation period T is (evaluating the integral in a 
saddle-point approximation) 



AE'"'{t) 



dt 



8a° 



bb rm^o 



h ujq + uj'^{4:XQmuJo/h ~ I) 



(22) 



If the particles are in different internal states, the center 
of mass does not decouple from the relative motion. No 
analytical solution is found in this case, and one must 
resort to numerical techniques to evaluate the coUisional 
phase 
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a) 





b) 




FIG. 2. Configuration at times t < 0, t > t (a) and at 
< t < T (b). The solid (dashed) curves show the potentials 
for particles in the internal state \a) {\b)). 



4-. Implementation of conditional dynamics 

If at time r the atoms have come back to their ini- 
tial spatial distribution, corresponding to a symmetrized 
product of the ground states of the two wells, then af- 
ter the barrier is raised they will remain trapped around 
the original position. The only change in the overall 
state will b e a phase (f)^^ + cj)^^ + as discussed in 

Sect. [IB 2 , Therefore, the gate operation time r has 
to be chosen in such a way as to maximize the over- 
lap |(VA'^=^(r)|V''^i''2(0))|2 for aU /3i,/32. If the modifica- 
tions in the atomic motion due to interaction are not too 
strong, this condition will be satisfied to a good approx- 
imation after an integer number n of oscillations. Thus, 
for r w nT , the following mapping is realized: 



\a)i\a)2 
\a)i\b)2 
ma)2 
mb)2 



Wi\b)2, 



(23) 



where = as discussed in Sect. |II C 2| . If we apply a 
further single-bit rotation |0)(0|e-*'^° + |l)(l|e~*('^'+'^°'') 
(where the logical states are defined as |0) = \a) and 
|1) = 1 6)) and take into account that for symmetry rea- 
sons cj)%^ = (j)^^ , the mapping Eq. (^3|) realizes the funda- 
mental phase gate 



|0)|0) 
|0)|1) 
|1)|0) 

11)11) 



|0)|0), 
|0)|1), 
|1>|0), 



-20" 



'|1>|1>, (24) 

where the phase difference (ff — 2(/)^^ has to be adjusted 
to ±7r by a proper choice of the trap parameters. 



III. PHYSICAL REALIZATIONS 

A physical implementation of the scenarios described 
m Sec.g requires an interaction which produces internal- 
state-dependent conservative trap potentials and the pos- 
sibility of manipulating these potentials independently. 
The choice of the internal atomic states \a) and \b) has 
to be such that they are elastic (i.e. the internal states do 
not change after the collision). To achieve entanglement 
operations with high fidelity, one has to be able to load 
or cool the atoms to the ground states of the trapping po- 
tentials. Finally, for the application of parallel quantum 
computing one needs periodic structures (e.g. optical lat- 
tices), together with the ability to control the positions 
of the atoms and to fill the lattice sites selectively. 
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A. Two-qubit gates in optical lattices 

In this Section we want to discuss how a number of 
difficulties can be overcome that one encounters when 
trying to use optical lattices for quantum computing. We 
will first show how one can achieve a filling factor of 1 
with particles in the ground states (lowest band) of the 
lattice. This can be achieved by using an ultracold very 
dense sample of weakly interacting atoms, namely a Bose- 
Einstein condensate, and slowly turning on an optical 
potential. The repulsive interaction between the particles 
increases as the optical potential is made deeper. At the 
same time the hopping rate at which particles move from 
one site to the next decreases. If the optical lattice is 
turned on on a time scale much slower than the hopping 
rate and the temperature kT can be kept much smaller 
than the interaction energy between two particles in one 
site, one can achieve a filling of the optical lattice with 
exactly one particle per lattice site, pq] Finally, we note 
that a filling factor of one out of two lattice sites has been 
achieved in very recent optical lattice experiments. [ |65| 

We will also discuss how the lattice potentials can be 
moved in a state-selective way for implementing the two- 
qubit gate . For alkali atoms with a nuclear spin equal 
to 3/2 we show how atoms in different hyperfine levels 
can be moved into different directions. It is clear that 
other difficulties like e.g. addressing single qubits exist, 
but they will not be discussed here since their experimen- 
tal solution is not specific to the present implementation. 



1. Hamiltonian for a Bose-Einstein condensate in an optical 

lattice 



We assume a Bose-Einstein condensate of atoms in in- 
ternal state I a) to be loaded into an optical lattice po- 
tential Vt{x) + Vb(x), where 

Vb(x) = sin2(fca;) -|- VyQ sin^ (ky) + V^q ain^ (kz) (25) 

is a periodic optical lattice potential and Vt(x) is a su- 
perlattice potential slowly varying in space compared to 
Vq (x) . k is the wave number of the lasers producing the 
lattice potential. The Hamiltonian reads |^ 

f d3x^t(x) (-^^^ + yo(x) + Vt{^) - V(x) 



d3a;?/;^(x)V'^(x)V'(x)V'(x), (26) 



where ipi^) is the bosonic field operator and fi is the 
chemical potential i.e. a Lagrangian multiplier to fix the 
number of particles. Expanding the field operators in the 
Wannier basis while keeping only the dominant terms fz^ j 
Eq. (pq) reduces to the Bose-Hubbard Hamiltonian 



H 



hi(hi 



<«j> 



-1), 
(27) 



where the operators hi = h\hi count the number of 
bosonic atoms at lattice site i; the annihilation and cre- 
ation operators hi and b\ obey the canonical commutation 



J is the tunneling matrix element 



relations fcj] = 
and U describes the (repulsive) interaction between par- 
ticles at the same lattice site, = Vxixi) is the value 
of the slowly varying superlattice potential at site i. The 
ratio U/J is controlled by the depth of the optical lat- 
tice potential Vjo . Increasing Vjo (via the intensity of the 
trapping lasers) reduces the tunneling matrix element J 
and increases the repulsive interaction between the atoms 

u M. 



2. Loading the lattice 

In order to perform gate operations in optical lattices 
we have to be able to selectively fill the lattice sites with 
exactly one particle. This can be achieved by making use 
of the phase transition from a superfluid BEC phase to 
a Mott insulator (MI) phase at low temperatures, which 
can be induced by increasing the ratio of the onsite in- 
teraction U to the tunneling matrix element J predicted 
by the Bose-Hubbard model jT^J?^. In the MI phase 
the density pi (occupation number per site) is pinned at 
integer n = 0, 1, 2, . . . corresponding to a commensurate 
filling of the lattice, and thus represents an optical crystal 
with diagonal long range order with period imposed by 
the laser light. Particle number fluctuations are thereby 
drastically reduced and thus the number of particles per 
lattice site is fixed. The number of particles per lattice 
site depends on the chemical potential in the isotropic 
case Ci = fl^ - In the non-isotropic case we may view 
fi — ei as a local chemical potential. Therefore pi can be 
controlled by the superlattice potential Vt(x). 



a) 




b) 
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FIG. 3. Superlattice potential in 2D 

with VT{x,y) = 40J (sin^(7ra;/9a) + sm^{TTy/9a)) with the 
spacing between the lattice sites. The particle density p{x,y) 
for four superlattice wells is shown. Parameters: a) [/ = 30J 
and p = 15J, h) U = 50J and p = 27 J. 



Using a Gutzwiller ansatz ||75|,[78|,^ for the wave func- 
tion we have performed a mean field calculation to 
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demonstrate how, by a proper choice of the potential 
Vt(x), one can fiU certain blocks of the optical lattice 
with exactly one particle at temperature T — 0. Figure 
^ shows the result of this mean field calculation, a MI 
phase where the lattice sites are either filled with or 1 
particles. The number fluctuations are almost equal to 
zero and thus not shown in this plot. To achieve a MI 
phase at finite temperature T 7^ one has to fulfill the 
requirement kT <^ U where the interaction strength U 
gives the order of magnitude of the first excitation energy 
in a MI phase. One also has to ensure that particles do 
not move from a filled site with energy to an adjacent 
empty site with energy Cj i.e. the temperature has to be 
much smaller than the energy difference between these 
two sites kT <^ ej — e^. In Sec. we will need periodic 
fillings of optical lattices as shown in Fig. ^ to implement 
efficient multi-particle entanglement operations and for 
parallel quantum computing. 



that the ac-Stark shifts of these two levels cancel each 
other. The optical potentials for these two states are 
K«.=±i/2 (z, 0) = a|^o P sin^ {kz ± 9) . 




FIG. 5. Level scheme of Rb and Na and laser configura- 
tion, (a) Fine structure energy levels and laser configuration. 
The detuning is chosen such that the polarizabilities a+_ and 
a — I- vanish, (b) Hyperfine level structure. 



3. Moving the lattice potentials state selectively 

We consider the example of alkali atoms with a nuclear 
spin equal to 3/2 (*^Rb, ^"^Na) trapped by standing waves 
in three dimensions and thus confined by a potential of 
the shape as given in Eq. (p5|). The internal states of 
interest are hyperfine levels corresponding to the ground 
state 5*1/2 as shown in Fig. ||b. Along the z axis, the 
standing waves are in the linZlin configuration (two lin- 
early polarized counter-propagating traveling waves with 



the electric fields Ei and E2 forming an angle 26 |80|) as 
shown in Fig. ^. 



26x 



20 

4^ 





El E2 

FIG. 4. Laser configuration along the z-axis. 



The total electric field is a superposition of right and 
left circularly polarized standing waves (cr^) which can 
be shifted with respect to each other by changing 9, 

E+{z, t) = Eqc-'-''^ [e+ sin(fcz+6') + el sin(fcz-6')] , (28) 

where e± denote unit right and left circular polariza- 
tion vectors, k — v/c is the laser wave vector and i?o 
the amplitude. The lasers are tuned between the P1/2 
and P3/2 levels so that the dynamical polarizabilities 
a±T of two fine structure S1/2 states corresponding 
to nis ~ ±1/2 due to the laser polarization cr^ vanish 

(a^ = a [- = 0), whereas the polarizabilities a±± due 

to (T^ are identical («++ — a_ = a). This configura- 
tion is shown in Fig. |^a and can be achieved by tuning 
the lasers between the P3/2 and Pi/ 2 hue state levels so 



We choose for the states \a) and \h) the hyperfine struc- 
ture states \a) = |F = 1, m/ = 1) and |6) = |F = 2, m/ = 
2) . Due to angular momentum conservation, these states 
are stable under collisions (for the dominant central elec- 
tronic interaction |^ , ^ ). The potentials "seen" by the 
atoms in these internal states are 

V{z, 9) = [Kn.=i/2(^, 0) + 3F„,=„i/2(z, 9)] /4 (29a) 
V\z,9) = V„,^^i,2{z,9). (29b) 

If one stores atoms in these potentials and they are deep 
enough, there is no tunneling to neighboring wells and 
we can approximate them by harmonic potentials. By 
varying the angle 9 from 7r/2 to 0, the potentials 
and V° move in opposite directions until they completely 
overlap. Then, going back to 9 = 7r/2 the potentials re- 
turn to their original positions. The shape of the po- 
tential V"" changes as it moves. By choosing 9{t) = 
^ (1 - (1 + eM-in/Tr?)) I (1 + exp((t2 - )/r2))) /2 
with = 25/a; and Ti = 25/aj, the frequencies and 
displacements of the harmonic potentials approximating 
are exactly those plotted in Fig. ^a. 



^. Gate fidelity 

We use the minimum fidelity F to characterize the 
quality of the gate. F is defined as 

F = min(^|trext ® Pc^tU^) (30) 

where \ip) is an arbitrary internal state of both atoms, 
is the state resulting from \(p) using the mapping (|l^). 
The trace is taken over motional states, U is the evolution 
operator for the internal states coupled to the external 
motion (including the collision), and pcxt is the density 
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operator corresponding to both atoms being at a tem- 
perature T at time t = -t In Fig. §3 the fidehty F 
is plotted as a function of the temperature T for the dis- 
placements and trap frequencies shown in Fig. This 
figure shows that one can achieve very high fidelities in 
realistic situations. 




FIG. 6. a) Upper plot: Displacements 5x°'{t)/d (solid line) 
and 1 + 5x''{t)/d (dashed line). Lower plot: Trap frequencies 
Lu''{t)/Lo (solid line) and u!^{t)/u! (dashed line), b) Fidelity F 
against temperature kT/huj for *^Rb with as — S.lnm. Here 
lu = 2tt X lOOkHz and d = 390nm. 



B. Two-qubit gates in magnetic microtraps 

We now consider the implementation of a switching po- 
tential by means of electromagnetic trapping forces. We 
first discuss the possibility of obtaining the desired state 
dependence by assuming some improvements on devices 
which are now experimentally available [^-^ . Then we 
compute the performance of a quantum gate for realistic 
trapping parameters. 



1. Microscopic electromagnetic trapping potential 

The interaction between the magnetic dipole moment 
of an atom in some hyperfine state \F,mp) and an ex- 
ternal static magnetic field B entails an energy t/magn ~ 
gpfJ-BTnElBl depending on the atomic internal state via 
the quantum number mp (here /is is the Bohr magne- 
ton and gp is the Lande factor). On the other hand, 
the Stark shift induced on an atom by an electric field E 
gives a state-independent energy Ud ~ ^a\E\'^, where a 
is the atomic polarizability. The interplay between these 
two effects can be exploited in order to obtain a trapping 
potential whose shape depends on the atomic internal 
state. As an example, we consider an atomic mirror with 
an external magnetic field ||8^|8^ , providing confinement 
along two directions with trapping frequencies which can 
range from a few tens of kHz up to some MHz. Micro- 
scopic electrodes can be plugged on the mirror's surface 
prt , thus allowing for the design of a potential with the 
characteristics described in Sect. H C 



2. Loading and moving atoms within the trap 

Several schemes of loading atoms into the trap have 
been envisaged (see for example [p4|-|86|). Most of them 
rely on an intermediate stage where atoms can be trapped 
and cooled without coming in contact with the magnetic 
mirror. This pre-loading trap can be either initially dis- 
placed from the surface, or close to it but based on a 
different trapping mechanism (for instance an evanes- 
cent wave mirror, where different internal states can be 
trapped by gravity before the atoms can be put in 
the right states for magnetic trapping) , to be replaced by 
the electromagnetic microtrap with a gradual switch-on 
of the electric and bias magnetic fields in the final stage 
of loading . This could also allow for implementing 
a controlled filling of the trap sites, in a similar way to 
that already discussed in Sec. 



HI A. A further feature 



to be implemented in view of performing more complex 
algorithms is the arrangement of several gate potentials 
in a periodic pattern, and the possibility of transport- 
ing atoms within this structure. An example would be 
given by two adjacent rows of potential minima, shiftable 
with respect to each other, where atoms co uld be lo aded. 
A system like the one suggested in Sect. IIIB 1 could 
allow in principle to obtain such a configuration, since 
the magnetic field minima can be shifted parallel to the 
surface by rotating the bias magnetic field. In this way 
it should be possible to move some atoms, while hold- 
ing others in place by means of additional local electric 
fields ]8^-|8^. Provided that atoms can be addressed in- 
dividually, which is needed even for performing a one-bit 
quantum gate, a procedure for implementing a simple 
quantum algorithm could be the following: perform a 
gate between two suitably chosen atoms, being close to 
each other but belonging to different rows, then mutu- 
ally displace the rows and select another pair of atoms, 
including one of those coming out from the previous gate. 
Repeat until the algorithm has been operated, applying 
the required one-qubit rotations in between the above 
steps and possibly performing some of them in parallel. 



3. Switching the trap potentials state selectively 

We choose for the states \a) and \b) the same hyper- 
fine structure states of ^^Rb considered in the previous 
Section, which are low-magnetic field seekers. If both 
particles are in state |a), there is no interaction-induced 

The re- 



phase shift, as already discussed in Sec. II C 2. 
suits for both particles in state |&) are shown in Fig. |^. 
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FIG. 7. Dynamics during gate operation, with both atoms 
in state \b): a) interaction- induced phase shift - the circles re- 
fer to the perturbative calculation (p^; b) projection of the 
evolved state on the corresponding state evolved without in- 
teraction; c) projection of the evolved state on the initial one. 
We choose ui ~ 23.4 kHz and uiy = lOz = 150 kHz, correspond- 



ing to ground-state widths ~ 50 nm, ay 



28 nm, 



with the initial wells having frequency loo ~ 2uj and displaced 
by xo = 3\/2ax. We take for the scattering length the known 
value for *^Rb, i.e. a'* ~ a^'' = 5.1 nm. Time is in units of 
the oscillation period T. 



The time dependence of (j)''^ is step- like (Fig. 0a): the 
collisional phase is incremented at times tk = (2fc— l)T/4, 
when the atoms meet at the center of the well, and re- 
mains constant at intermediate times, when they sepa- 
rate again. The influence of the interaction on the atomic 
motion can be seen from Fig. depicting the overlap 
between the evolved interacting two-atom state \^''^{t)) 
and the corresponding state |V'(o)(0) computed without 
taking into account the interaction. The curve has local 
minima at times tk, signalling that a collision is taking 
place, and shows a global decrease corresponding to a 
slight delay of the interacting motion with respect to the 
non-interacting one. As it can be seen from Fig. |7|3, this 
effect is not dramatic: the oscillation period in the pres- 
ence of interaction is increased just by ST « 2 x 10~^T 
(with the parameters used here), and the harmonic po- 
tential ensures that the system comes periodically back 
to its initial state. After 7 oscillations we get a phase 
shift due to the interaction of tt, whereas the pertur- 
bative formula (^2|) gives 7(f>^ « 0.987r. Therefore we 
choose r = 7(T -t- ST) w 0.15ms: the overlap between 
the initial and the evolved wave function at that time 
is |(V'''''(t)|V'''''(0))P « 0.996. The behavior turns out to 
be quite different when the atoms are in different 
internal states: the phase shift increases more rapidly, 
but after a few oscillations the system does no longer 
come back to the initial state. This has a simple expla- 
nation. The two atoms collide as soon as the one being in 



state \b), moving within the potential mw^a;^/2, reaches 
its turning point, where the other atom is trapped. The 
interaction time is therefore longer than if both atoms 
were in state \b). Indeed, in that case they meet at the 
trap center, with their maximal velocity. This explains 
why the system picks up a bigger phase shift per oscil- 
lation period in the present case. On the other hand, 
the collision excites the motion of the atom in state \a) 
within its own well, and therefore the initial state is no 
longer recovered. This problem can be avoided if the 
potential minimum for state \a) is displaced along the 
transverse direction from the one for state \b) by means 
of an additional electrostatic field p^-|8^, so that the 
atoms interact if and only if they are both in state \b). 
This problem would not exist in an adiabatic scheme for 
the gate operation, when the shape of the potential is 
changed slowly with respect to the atomic motion. This 
will be the subject of future investigation. 



4. Gate fidelity 

The calculation of the fidelity in this case has to take 
into account the symmetrization of the wave function un- 
der particle interchange, expressed by an operator S to 
be explicitly inserted into Eq. (^0[): 

F = mm {tre,t i\ip) {ip\ ^ Pext) SW\^)] } , (31) 

With the parameters quoted above, we obtain F > 0.98. 
In order to reach such a fidelity, timing has to be quite 
precise, with a resolution of the order of 10~^T corre- 
sponding to tens of ns in this case. 



IV. PARALLEL QUANTUM COMPUTING 

In this Section, we will discuss how quantum gates 
based on controlled collisions can be exploited for quan- 
tum computing. It is clear that, with the realization of a 
universal two-bit gate, any quantum computation can be 
performed, just as it is the case with other implementa- 
tions. On the other hand, manipulations such as moving 
and switching potentials offer a great deal of parallelism 
p3| , p[ not available in other systems. 

We will focus our attention on implementations in opti- 
cal lattices. Some of the ideas could readily be translated 
into arrays of magnetic microtraps, if the distances be- 
tween the individual potential wells could be made much 
shorter than present-day state-of-the-art of nanofabrica- 
tion. In such a situation, adiabatic variants of the switch- 
ing operations (see comment at the beginning of Sec. |II C| ) 
can be used to create multi-particle entangled clusters 
of neighboring atoms, similar as with moving potentials. 
Details of this analysis will be presented somewhere else 
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One may ask, what can be done in optical lattices that 
cannot be done in other implementations? The answer 
to this question depends on a number of experimental 
conditions such as the possibility of creating regular fill- 
ing structures and, like in ion-traps, on the possibility of 
addressing sin gle at oms individually. In the following, we 
will first (Sec. [V A) give an example of what can be done 
with controlled lattice movements in conventional set-ups 
i.e. with random filling of the lattice sites and without 
any control of the position of individual atoms. We will 
see that this already allows one to perform interesting 
spectroscopic studies of the degree of ent angle ment be- 
tween the atoms thus created. Next (Sec. IVB), we will 
describe what can be done if one achieves a regular occu- 
pation of the lattice sites and can address the atoms in- 
dividually. Under such circumstances, an efficient imple- 
mentation of quantum error correction and of a quantum 
memory (concatenated Shor code) is possible. Further- 
more, fault tolerant versions of certain quantum gates 
and of quantum error correction can be im plem ented 
straightforwardly, as will be sketched in (Sec. IV C). Fi- 
nally, in Sec. IV D, we describe how auxiliary atomic lev- 



els can be used to realize highly selective entanglement 
operations, where individually selected atoms are swept 
across the lattice to create GHZ states [M of a large 



number of particles. Together with IVB and IV C , this 



scheme has all the ingredients that are necessary for an 
efficient realization of fault-tolerant quantum computing. 



A. Multi-particle entanglement operations 



The two-qubit gates described in Sec. Ill correspond 
abstractly to an atom interferometer as shown in Fig. ||. 
The interferometer has two inputs which are the two 
atoms trapped at neighboring potential wells. By shifting 
the potentials back and forth as described in Sees. II B, 



only one combination of paths of the two particles over- 
laps and leads to a phase shift, namely the paths corre- 
sponding to state |a)i for the left particle and \h)2 for 
the right particle. To emphasize the role of the internal 
states as logical states, we shall henceforth use the nota- 
tion |0) = \a) and |1) = \b) and neglect the kinetic phases 
(ff , (ffi as they appear in (p^). Furthermore, we drop the 
atomic index as long as there is no danger of confusion. 



(a) 



time 




(b) 



out 



10)10) 10)10) 

10) 11) 10)11) e'* 

11) 10) 11)10) 

11)11) 11)11) 

atom 1 atom 2 
FIG. 8. Atom-interferometric process realizing the quan- 
tum gate, (a) Two-particle interferometer; (b) Truth table. 



The logical truth table corresponding to the interfero- 
metric process is shown in Fig. [A similar identifica- 
tion of logical states can be made in magnetic traps as 



is pointed out in Sec. II C 4. The labelling of the paths 



for the left particle in the interferometer has to be inter- 
changed in this case.] For — cffi^ = tt this realizes a 
phase gate ^ . The phase gate and the set of all one-bit 
unitary transformations, which can be realized by Ra- 
man laser pulses on the internal states 10) and |1), define 
a universal set of quantum gates. pl|-p^ 

An important difference between optical lattices and 
other implementations is given by the global effect of the 
lattice manipulations. To illustrate this point, consider 
first a two dimensional lattice as in Fig. ^ with random 
occupation of the sites and a filling factor 77 ^ 1, where 
Tj is defined as the average number of atoms per lattice 
site. Let us assume that the loading of the lattice can be 
accomplished in such a way that there are no multiply 
occupied lattice sites, i.e. that each lattice site is occupied 
by no more than a single atom. 

Selected lattice region (2-dim): 





triplets 



pairs 



isolated 
atoms ~ r| 

FIG. 9. Random occupation of a two dimensional lattice 
with single atoms. 

Then, in any region of the lattice, one will find isolated 
atoms, pairs of neighboring atoms, triplets, and so forth, 
with a relative frequency proportional to rj, , rf , respec- 
tively. Consider now the following Ramsey experiment 
psf where initially all atoms are prepared in the inter- 
nal state |0) and in the motional ground state of their 
individual potential wells. In some selected region of the 
lattice, the following sequence of operations is applied: 
(1) a 7r/2 laser pulse brings all atoms into a superposi- 
tion of the internal states |0) and |1); (2) the lattice is 
shifted across one lattice site and then, after a variable 
length of time, shifted back to its original position, (3) 
finally a second 7r/2 pulse is applied to the region. The 
effect of this sequence is illustrated in Fig. |l^. For a 
group of iV — 1,2,3, .. . neighboring atoms, the lattice 
shift corresponds to a N-particle interferometric process. 
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1 e'* 

|BELL> 



1 + e'* 1 e'* 

\W) = 1000) + |GHZ) 

Et/H E/U Er/U 




|0> |0> |0> |0> |0> 

FIG. 10. Entanglement of pairs (left) and triplets (right) 
of neighboring atoms by a single lattice shift. 

Specifically, one obtains the following transformations. 
For isolated atoms: 



|0) — |0); 
for pairs of neighboring atoms: 



100) 



1 



-100) 



1 



ai4> 



2 ' ' 2 
and for triplets of neighboring atoms: 



1 000) 



1 



-|000) 



1 



2 ' ' 2 
where we have used the notation 
1 



I bell) ; 



GHZ) 



(32) 



(33) 



(34) 



BELL 



V2 
1 



{|o>|+>-|i>|->}, 



|GHZ) = ^{|0)|+)|1)-|1)|-)|0)}, 



(35) 



and |±) = (|0) ± |l))/-\/2. The expressions for groups 
of more particles become more complicated and shall be 
ignored in the present discussion. It is clear that for 
(p = TT Bell- and GHZ states p6| , po| are created by a sin- 
gle lattice shift at various places within the region. This 
corresponds to an ensemble of 2-bit and 3-bit quantum 
gates, respectively, acting simultaneously at different lat- 
tice sites. To analyze the states ( ^ and (34) spectro- 
scopically one could measure the state of the atoms in a 
final step of the above Ramsey sequence e.g. by a fluo- 
rescence measurement. It is clear that by such a mea- 
surement the entangled states will be destroyed. On the 
other hand, by repeating this sequence many times with 
different samples, one can measure the fluorescence sig- 
nal as function of the phase cj) (interaction time). Under 
ideal circumstances, all isolated atoms will remain in the 
dark state |0) while all fluorescence signals come from 
Bell (~ 77^) or GHZ (~ 77^) states [0. To check that en- 
tangled states, rather than mixtures, are created in the 
process, the experiment is performed with different in- 
teraction times, e.g. times corresponding to (/) — tt and 
(j) — 2tt. For entangled states as in ( |3^ ) and ( |3^) all 
ffuorescence signals will vanish at (/) = 27r, while this will 
not be the case if the states created by the atomic colli- 
sions are mixtures of classical many-particle states. More 



generally, by measuring the visibility of the fluorescence 
signal one may study the fidelity of the entanglement cre- 
ated in the process, and its dependence on certain noise 
sources such as a finite temperature of the atoms. This 
way, the curve plotted in Fig. ||b) could be tested exper- 
imentally. 



B. Quantum error correction 

To employ these entanglement operations for quantum 
computing, one has to have precise control over the num- 
ber and the location of atoms that are involved in the col- 
lisional process. In addition to the ability of addressing 
single atoms, one therefore has to achieve a certain or- 
dered occupation of the lattice sites. As described in Sec. 
HI. A., Fig. ^, this can be by achieved by controlling the 
intensity of the trapping laser at sufficiently low temper- 
atures. This way optical crystals with periodic patterns 
of atoms can be created as indicated in Figs. |ri| and |l^ 

M 



FIG. 11. Ordered arrangement of atoms in an optical lat- 
tice (see also Fig. 



Under such circumstances the parallelism of the lattice 
manipulations can be exploited advantageously. On one 
side, similar logical operations can be performed simul- 
taneously at different locations on the lattice. On the 
other side, as we have seen in Fig. |l^, a single lattice 
shift can entangle whole groups of atoms. Two typ es of 
such entanglement operations are shown in Fig. |lj. One 
involves only the logical states |0) and |1), while the sec- 
ond uses a third atomic level as a "transport state" (see 
Sec. IVD), into which any atom must first be activated, 
before it can participate in an entanglement operation. 
In the following, we will first discuss applicatio ns of the 

1(a). 



IV D 



shift operation as in Fig. |lj(a). Later, in Sec. 
will consider a more flexible ("sweep") operation shown 
in Fig. |l|(b). 
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(a) 




(b) 





|r> 
transport 
state 



ID- 
FIG. 12. (a) "Shift operation": The internal atomic states 
|0) and 1 1) couple to different lattice potent ials tha t are moved 
against each other as explained in Sec. 








[II A 3 



This cor- 
responds to a multi-particle interferometer where the same 
phase shift is acquired whenever two paths temporary over- 
lap. By simple lattice manipulations, therefore, entire groups 
of atoms become entangled, (b) "Sweep operation" : For more 
selective entanglement operations, a third atomic level |r) is 
used In this scheme, only atoms in the level |r) are 

moved, whereas the states |0) and |1) are kept in the same 
potential. At the beginning of an entanglement operation, 
the atoms are first excited from one of the states |0) or |1) 
to the state \r) before the lattice is moved. This scheme is 
much more selective in the sense that those atoms which shall 
participate in a gate operation are first activated, before they 
couple to the moving lattice, and the collisional phases <j!>j can 
be varied for each interaction individually. 



consist of products of certain GHZ states. Abstractly 
speaking, the encoding operation consists of a mapping 
(embedding) of the qubit's 2-dimensional Hilbert space 
H into a 2^-dimensional Hilbert space of the form 



n 9 a|0) + /3|1) ^ a|Os) + /3|ls) e He C K 



(37) 



To stabilize the encoded information against decoher- 
ence, the code must be measured and corrected on a time 
scale r ^ 1/87 where 7 is the rate of decoherence for a 
single qubit. This is possible since all errors that may oc- 
cur on any one of the qubits of the codewords ( ^ ) map 
the code into a family of 2-dimcnsional subspaces of 7i®^ 
which are all orthogonal on He @1 ■ 

The Shor code (|36|) can be implemented efficiently in a 



two dimensional lattice configuration! IOC] as in Fig. |T|, 
by using the shift operation of Fig. |12|(a). To see this, 
imagine that the qubit/atom whose state is to be encoded 
is surrounded by neighboring atoms as in Fig. |l^. The 
idea is of course to encode the central qubit in the whole 
block of 3 X 3 qubits. Initially the central atom is in the 
unknown state jV') — ct\0) + f3\l) while all neighboring 
atoms are in state |0). As is shown in Appendix B, the 
initial state is transformed into the Shor code by a simple 
sequence of horizontal and vertical lattice shifts combined 
with certain 1-bit rotations, as indicated in Fig. ^(a). 
By this process, the information contained in "0 is so to 
speak de-localized over the whole block of 9 atoms. To 
check whether an error has occurred on one of the qubits, 
the block is first decoded by the inverse transformation 
pO[ , which involves the same sequence of lattice shifts as 
the encoding. Subsequently, one measures which of the 
neighboring atoms are in the state |1). 



(a) 



Encode 



V = alO> +|3I1> 



aio:ij+pii:^ 

Shor Code 



An application of shift operations as described in 
Fig. ^(a) concerns the reahzation of a quantum memory, 
where a qubit a\0) + /3|1) (with unknown coefficients a 
and f3) is encoded in the quantum state of a larger block 
of atoms and stabilized against decoherence with the help 
of quantum error correction p7|-p^. A particular quan- 
tum code that is able to protect a qubit against general 
1-bit errors (spin flip and phase flip) has been proposed 
by Shor It is a 9-bit code where the codewords 

|0s) = 2-3/2(1000) -I- 1111))(|000) + |111))(|000) + |111)) 
lis) - 2-3/2(1000) - |111))(|000) - |111))(|000) - |111)) 

(36) 



li>=^r-T--.. ^ " ^ . / . 

• [ ' ' Decode 

V^alO)'-pil) 

FIG. 13. (a) Encoding of a qubit into a block of 3 x 3 atoms; 
(b) Decoding and syndrome measurement. 

In the language of quantum error correction, the sur- 
rounding atoms of the central qubit in Fig. |l3|(b) are the 
carriers of the error syndrome |50[| , meaning that their 
state gives information regarding what type of error oc- 
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curred and, more importantly, which unitary 1-bit rota- 
tion has to be applied to the central qubit to restore it to 
the original state. In a fluorescence measurement, this in- 
formation corresponds to a specific pattern of bright and 
dark atoms surrounding the central qubit. For example, 
in Fig. |l3|(b) a spin flip has occurred in the central qubit. 
This means that the state of the block is transformed 
into Q!O-^'''|0s) -1-/3(7^''' I Is) where crj;''' is the corresponding 
Pauli spin operator associated with the central atom. By 
the decoding operation of Fig. ^(b), this state is mapped 
into a product state of 9 qubits in which the neighboring 
atoms in the central row are both in the (bright) state 
|1), while the other surrounding atoms are in the (dark) 
state |0). The central atom is in a state equivalent to 
\ip) up to a unitary transformation. Similar fiuorescence 
patterns are obtained if a phase error occurs in the cen- 
tral atom or an arbitrary 1-bit error on any of the atoms 
in the block. A complete table of the error syndrome is 
given in the Appendix. 

The essential point is that the measurement on the 
surrounding atoms does not reveal nor destroy the state 
of the central atoms (the coefficients a and P remain 
unknown throughout the process). If the sequence of op- 
erations "decode-correct-encode" is repeated sufficiently 
often within the decoherence time I/97 of the block, the 
state ip may be protected over arbitrarily long times, in 
principle. Here one assumes, of course, that the decod- 
ing and encoding operations themselves are free of er- 
rors. In our situation this means that all phases ac- 
quired in the atomic collisions can be perfectly controlled. 
Since these operations will always bear some imperfec- 
tion/imprecision, the probability that an error is intro- 
duced by an imperfect operation increases/accumulates 
with repeated applications of these operations. The gen- 
eral solution to this quantum-memory problem was given 
by Knill and Laflamme and by others [5^-|6^, and 
requires a concatenation of encoding operations as shown 



schematically in Fig. 14 




FIG. 14. Concatenated quantum coding. At each coding 
level, a single qubit is encoded in a block of a larger number 
of qubits, here 5. (See e.g. [p2[). 

The number of required concatenation steps depends 
on how long the qubit is to be stored. It can be shown 
that, given the precision of the operations is above a 
certain threshold, a qubit can be stored for an arbitrary 
long time, where the number of qubits required for en- 



coding (i.e. the length of the code) grows polynomially 
with the length of the storage time. In the optical lat- 
tice configuration, a concatenation of the encoding can 
be implemented straightforwardly. Imagine that, in the 
central block in Fig. O, the center atom is initially in 
state \ip) ~ a\0) -f /3|1) (similar as in Fig. |l^) whereas all 
other atoms are in state |0). This means that both the 
surrounding atoms in the center block and the atoms of 
all the other blocks are initially in state |0). The first 
step of the encoding operation is identical as in Fig. |l^ 
and results in the configuration where the center block 
is in a superposition of the Shor code words |0s) and 
I Is), whereas the surrounding blocks remain in state |0). 
In the second step, the same operation is repeated on 
a larger scale, i.e. the lattice is shifted across a larger 
distance such as to make the blocks temporarily overlap 
while the 1-bit operations of the first step are now re- 
peated on corresponding atoms of the outer blocks. As a 
result, the information \ip) originally carried by the cen- 
ter atom is now delocalized over 9x9 = 81 atoms! This 
scheme may be iterated as indicated in Fig. n^. 



FIG. 15. Concatenated quantum coding in an optical lat- 
tice. At each coding level, two-dimensional lattice displace- 
ments with an increasing periodicity are applied. The nested 
character of Fig. |l^ is here reflected by a self-similar filling 
pattern of the lattice. 



When in the second (and higher-order) encoding step 
the blocks are brought to overlap, one has to make sure 
that only phases between corresponding atoms of the dif- 
ferent blocks are accumulated. The most elegant way 
to achieve this would be with the aid of a technique 
where the and 1 states are displaced vertically before 
the atoms are moved. This could be i mplem ented in a 
three-dimensional lattice configuration |101|. The shift 
operation is then really a "lift & shift" operation. The 
coUisional interaction is then only switched on by vary- 
ing the vertical displacement, after the blocks have been 
moved horizontally. If such a lifting technique can not 
be implemented, e.g. in a truly two-dimensional config- 
uration, then during the horizontal motion there will be 
also collisions between non-corresponding atoms, for ex- 
ample the atoms in the right column of one of the blocks 
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with atoms in the left column of a neighboring block. To 
avoid these unwanted phase shifts, it is possible to vary 
the velocity of the lattice movement in such a way that 
during unwanted collisions a phase of e^'^' is acquired. 
This method is clearly more susceptible to decoherence. 
On the other hand, our numerical studies have shown 
pTf , that by an appropriate choice of the displacement 
function 9{t) in Fig. ^, the phase of a single collision 
can, in principle, be controlled with a very high precision 
(with fidelity > 0.9997) and the probability for exciting 
phonons remains correspondingly small |102| ] . 

It does not seem impossible that 9{t) could be con- 
trolled precisely enough to meet the threshold of fault- 
tolerant computation [ ]62| , but we have not yet made 
detailed numerical investigations for this situation. In 
summary, the method of concatenated coding can be im- 
plemented in optical lattices by repeated sequences of 
lattice displacements on self-similar filling structures. 



the NOT gate, while for example the Hadamard trans- 
form is more involved and requires a measurement with 
auxiliary qubits. Whether or not one can find similarly 
efficient implementations for a complete set of fault tol- 
erant gates, is still under investigation. 

To check whether an error has occurred during a gate 
operation, one has to measure whether the blocks are 
still in a superposition of the correct codewords. For the 
Shor code, this can be done in the following way p7[ , 
see Fig. ^ To detect a spin-flip, one has to measure the 
parity of the first two atoms in any row and compare it to 
the parity of the last two atoms of every row. To do this 
one would use an "Armada" of 3 x 2 ancillas in the state 
(|00) -I- |11))(|00) -I- |11))(|00) +J11)), which approaches 
the block from the left in Fig. |l^ by moving the lattice 
horizontally. 



C. Fault-tolerant computing 

In a quantum computer, we do not only wish to store 
quantum information, but also to process it in a quan- 
tum algorithm. To prevent an accumulation of errors 
during the calculation due to imperfect gate operations, 
one needs to use fault-tolerant quantum gates that act on 
the encoded information. Furthermore, errors should be 
corrected fault tolerantly, that is, without decoding the 
information (and therefore exposing the qubit to decoher- 
ence). The general theory of fault-tolerant computation 
has been developed by several researchers |62). In opti- 
cal lattices, many of such fault-tolerant operations have 
a geometrically intuitive implementation. For example, 
if two qubits are encoded in blocks of 9 atoms each, as in 
Fig. |l3|, a controUed-NOT operation can be implemented 
by moving one block on top of the other so that each pair 
of corresponding atoms from the two blocks share a sin- 
gle potential well and acquire a phase shift e*'^ . [This is a 
straightforward generalization of the situation in Fig. p|. 



(a) 



(b) |Q)|Q') 


IQ>IQ'> 






IOs>lls> 




iis)ia) 









m iQ> 

FIG. 16. Implementation of a fault-tolerant CNOT gate. 

When a 7r/2 pulse is applied on one of the blocks before 
and after the blocks are shifted, a fault-tolerant realiza- 
tion of the CNOT gate, with a truth table as in Fig. |l^ 
is realized. The minus sign may be eliminated by apply- 
ing a 37r/2 pulse instead of the second 7r/2 pulse. Sim- 
ilarly, one can find a simple fault-tolerant realization of 



FIG. 17. Implementation of fault-tolerant error correction. 



To measure the parities, the Armada is moved on top 
of the first two columns of the data block so that the 
atoms interact pairwise with atoms of the data block 
and acquire a phase shift of e^^ . To satisfy the cri- 
teria for fault tolerance, we need to avoid collisions 
while the ancillas are moved on top of the code, and 
thus need a "lift & shift" implementation of the oper- 
ation, as mentioned earlier. Suppose there was a spin- 
flip in one of the atoms of the first row. Then the 
state of the ancillary atoms after the interaction reads 
(-|00)-|-|ll))(|00)-f|ll))(|00)-h|ll)), and the error will be 
detected by measuring the parity of the ancillas in each 
row, after applying a Hadamard transform. In a second 
run, the Armada is reset in the initial state and then is 
moved on top of the last two rows of the block, and so on. 
To detect a phase-flip, a similar procedure is used with 
an Armada of 2 x 3 atoms that approaches the block in 
Fig. ^ from below by moving the lattice vertically. Since 
these ancillas should measure any change of sign in any of 
the GHZ states that make up the codewords (36), they 
have to be prepared in the state 1 000000) -f|llllll).A 
phase flip can then be detected as previously, where now 
a Hadamard transform has to be applied to the block 
first, before the "attack" starts from below. In the spe- 
cific implementation using optical lattices, one could also 
think about other schemes using only a single row of an- 
cilla atoms on each side of the data block in Fig. |l^ as 
realized in Fig. ||b). 
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D. Selectivity and "sweep operations" 

The examples discussed so far make use of the par- 
allchsm of the lattice shift to implement certain multi- 
particle entanglement (or gate) operations efficiently. 
On the other hand, the shift operation as described in 
Fig. 12(a) is too rigid, when certain operations should 
apply to a selected group of atoms only. This problem 
can in principle be solved by using a third atomic level 
|r) as indicated in Fig.[l^(b). In this scheme, the level \r) 
couples dominantly to a transport lattice ||99|] , while the 
"logical states" |0) and |1) are kept in the same poten- 
tial. At the beginning of an entanglement operation, the 
atoms are first excited from one of the states |0) or |1) to 
the state |r), before the lattice is moved. This scheme is 
much more selective in the sense that those atoms which 
shall participate in a gate operation are first activated, 
before they can participate in the lattice movement. All 
operations that we have discussed can then be realized in 
the same manner, with the additional property that only 
those atoms, to which the operation |1) |r) is applied, 
will participate. With this additional feature, it is clear, 
that universal computations can be implemented. 



li>- 



sweep 



{IO>+|r» (|0)+I1»(I0>+I1» ™ (I0>+I1» 

FIG. 18. Realization of an (A''+ l)-dimensional GHZ state 
by a single sweep operation. 

Another merit of this scheme is that one can realize 
more flexible entanglement operations. Consider, for ex- 
ample, a 1-dimensional situation as in Fig. |l^ with a 
string of N atoms initially prepared in the product state 
(|0) + |1))®^ and a selected additional atom (left) in 
the state (|0) -|- |r)). By moving the transport lattice, 
the selected atom is swept across the N lattice sites. 
During that motion, it interacts with each of the N 
atoms thereby transforming the state of each atom into 
e^-^^lO) + e*"^']!), with a differential phase cj) = cj)^ - 0°. 
The resulting total state is of the form 



AN4,'- 



(|0) 



|0)(|0) + |1))(|0) + |1))...(|0) + |1)) 
e'^|l))(|0)+e'^^|l))...(|0)+e^^|l)). (38) 



As long as the coUisional phases are different {(p ^ 0) for 
the two logical states, (j) can by varied with the speed of 
this sweep operation. For (p = tt one obtains a -f 1 
dimensional GHZ state (see Fig. |l^ ) which can easily be 
brought to the standard form 

|^) = i=(|0)|0)|0)...|0) + |l)|l)|l)...|l)). (39) 

Note that for the creation of this state only a single sweep 
operation is required! 



This scheme can be generalized in several directions. 
By varying the speed by which the lattice is moved dur- 
ing the sweep operation, the phases can be controlled 
individually for each atom of the string as indicated in 
Fig. allowing for more complex entanglement oper- 
ations. As a final example consider a configuration as 
in Fig. |9(b), with a "source register" consisting of a 



string of m atoms in the state \a) = \ai a2 0,3 ■ ■ ■ Om), 
ttj £ {0,1} and a "target register" of m further atoms in 
the state (|0) -t- ll))*"", similar as in Fig. M 



sweep 



\a,a,a, .... «,„.,«„,) (IO>+|1»(IO>+|1» (|0>+ll» 

FIG. 19. Implementation of the quantum Fourier trans- 
form by a sweep operation with variable speed. 

The state vector |a) = |ai 02 03 • • • am) should be 
interpreted as a binary representation of the number a = 
ai2""i-|-022'""2-|-. . .-|-a„2°. Consider now the following 
operation where the source register is first activated to 
couple to the transport lattice, meaning that each of the 
atoms 1 to m that is in state |1) is excited to state |r). 
Next, the lattice is moved to the right so that atoms of 
the source and the target register interact; this motion 
continues with variable speed until the source register 
completely overlaps with the target register. It is helpful 
to mentally decompose this operation into discrete steps. 
In the first step, the transport lattice is shifted one lattice 
site to the right such that the mth atom of the source 
register interacts with the first atom of the target register. 
One can tune the interaction time such that a certain 
phase shift is acquired during this interaction, namely 
(j> — 271/2"^. In the next step, the transport lattice is 
moved one lattice site further to the right such that now 
the m th atom of the source register interacts with the 
second atom of the target register, while at the same time 
the m — 1 th atom of the source register interacts with 
the first atom of the target register. In this step, the 
interaction time is made double as long as in the first 
step, so that cj) = 27r/2™~^, and so on. After the lattice 
has been moved across m sites in this vein, the total state 
of the source and the target register is given by 

^2iTi0.aia2...ajj 



\ai 02 03 



®(|0) 



(|0) 



^27ri0.a2 ■ ..a 



(|o) 



1)) 



(40) 



Finally, the lattice is shifted back to the original posi- 
tion without changing the phases any more (modulo 2tt, 
see earlier remark, or the process can be made symmet- 
ric such that only half the phase values are accumulated 
during the motion to the right while the second halves 
of the phase values are accumulated when the lattice is 
brought back to its original position.) The overall effect 
of this sweep operation can be summarized in the form 

\n\\ \ > 



(41) 
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wherein \a) and | ) denote the initial state of the source 
and the target register, and \ J-{a)) is the discrete Fourier 
transform of the function j aj . The additional phase 
factor e**^°^ accounts for a possible phase shift arising 
from collisions among different atoms of the source, if 
no vertical displacement of the transport lattice is pos- 
sible. It should be remarked that for a superposition of 
different input states the source and the target register 
become entangled, and therefore a direct a pplication of 
this method in the Shor algorithm [ 103 104 1 is not possi- 
ble. Nevertheless, this final example demonstrates a re- 
markable flexibility of the entanglement operations that 
are possible in optical lattices and similar systems, of- 
fering new perspectives for efficient implementations of 
quantum algorithms. 



V. FINAL REMARKS 



also thank David DiVincenzo and Andrew Steane for 
helpful discussions on fault-tolerant quantum comput- 
ing during the Benasque Workshop 1998. One of us 
(T. C.) thanks M. Traini and S. Stringari for the 
kind hospitality at the Physics Department of Trento 
University, and the ECT* for partial support dur- 
ing the completion of this work. This work was 
supported in part by the Osterreichischer Fonds zur 
Forderung der wissenschaftlichen Forschung, the Euro- 
pean Community under the TMR network ERB-FMRX- 
CT96-0087, the Institute for Quantum Information 
GmbH, and by the Schwerpunktsprogramm "Quanten- 
Informationsverarbeitung" der Deutschen Forschungsge- 
meinschaft. 



APPENDIX A: ONE PARTICLE IN A MOVING 
HARMONIC POTENTIAL 



It is clear that, at the present time, most of the exper- 
imental requirements have yet to be realized, before one 
can implement quantum computing. There are, however, 
recent achievements in cooling and trapping of atoms in 
optical lattices and in magnetic microtraps which make 
it seem possible that some of these elements could be im- 
plemented in the laboratory in the near future. There are 
short-term and long-term perspectives. Essential for all 
quantum information experiments is a successful cooling 
of the atoms to the ground state of a three dimensional 
lattice. Numerical calculations |3^ using realistic param- 
eters give kT < 0.2fLU! as a critical value. Under these 
circumstances, one could perform interesting Ramsey- 
type spectroscopic studies of the fidelity of multi-particle 
entanglement as discussed earlier. To do this, neither 
single-atom addressability is required nor are regular fill- 
ing structures. When the latter requirements can be re- 
alized, on the other hand, coding experiments can be 
done and a quantum memory be implemented. Finally, 
if one can find three-level schemes with different scatter- 
ing phases for the logical states, universal computations 
can be performed. The parallelism of the lattice could 
then be exploited for efficient implementations of fault- 
tolerant quantum computing. 

We have discussed multi-particle entanglement 
schemes mainly in the context of optical lattice imple- 
mentations. Some of these ideas could readily be adopted 
in implementations with magnetic microtraps if one uses 
adiabatic schemes. A basic requirement for this is the 
possibility of creating quantum dots that are spatially 
sufficiently close to each other. These ideas will be dis- 
cussed somewhere else 1891. 
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1. Hamiltonian 

The center of the potential with frequencies uJx, and 
LOz is assumed to be given by x(t) — (a;(t),0,0) and the 
Hamiltonian reads 

H^H^ + Hy + Hz, (Al) 

where iJ^ = huJz{ctlo,z + ^/'^), Hy = hcuy {aj^ay + 1/2), and 

= hw., ^ala., + - + « + a,)-^ + j . (A2) 

The a's are bosonic destruction operators and x{t) is 
given in harmonic oscillator units. We will concentrate 
on the x-direction leave out the subscript x and normal- 
ize energies to huj. 

2. Exact solution 



The Schrodinger equation for the Hamiltonian Eq. A2 
can be solved exactly [105,73|. To do so wc define 



Hq — a? a + -, 



1 ft 



(A3) 



(A4) 



and 



/3{t,-T)=t f ds (K{s-T){dsK*{s-T)}+'-^ 



(A5) 

If initially at time — r the system is in the state 
I ^ ("''')) = where \n) is the n-th harmonic oscilla- 
tor eigenstate we get 
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The kinetic phase <j) is thus given by the phase of the 
overlap of with the instantaneous ground state 

£'(x(i))|0), where D{-^) — exp(7a^ — 7*0) denotes the 
displacement operator 



-argi 



|Z?(x(i))t|vI/(i))) 



(A7) 



The interaction phase can be found by Eq. (Q) with the 
known \^{t)). 



3. Corrections to the adiabatic approximation 

We assume x{t) to be an analytic function of t and 
that x{t) > dtx{t) > d'^x{t) > ... > d^x{t). By ex- 
panding in orders of the time derivatives we can write 
for K{t', -t) 

z"+i{a,"5(s)}e'(''+^)|^z*„^ , (A8) 

n=0 / 

where iV is a positive integer. Note that if we may ne- 
glect all terms of order greater than dtx{t) and start in a 
coherent state |*(-t)) = D[x[-T)+idtx{t)\t=^r)\^) the 
state will always be a coherent state with {x{t)) = x{t) 
and {p{t)) = dtx{t). 

Now we assume for simplicity that x{t) = x{—t) — 0, 
{dtx{t))\t=-T - {dtx{t))\t=r = and {d^x{t))\t^^r = 
(— l)"((9"a;(i))|(=7. for n > 1. The system is assumed to 
be in the state |^(— t)) = |0), initially. We keep all the 
terms to fourth order in the derivatives (in the integrand) 
and find 

K{t,-T) = i= [{^{^^, - dt)x{t)] (e'^'+^^ - 1)- 



and 



V2 

{93s(t)}(e^(*+-) + i; 



(A9) 



+ \{dlx{t)Y{l-e-'^'^^^ 



(AlO) 



If we choose (t + t) = 2n7r with integer n the largest 
correction to t he ap proximation to the kinetic phase dis- 
cussed in Sec. [IB 1 is of third order. Also the amplitude 



of the first excited state is of third order as can be seen 



from Eq. (A9) 



APPENDIX B: QUANTUM ERROR 
CORRECTION AND THE IMPLEMENTATION 
OF SHOR'S CODE 

Consider a one-dimensional configuration with a string 
of n atoms, where XQ.x1.x2, . ■ ■ ,Xn G {0,1} label the 
internal state of the atoms at position 0, 1, 2, . . . , n of the 
lattice. An elementary lattice-shift operation as given in 
Fig. |lj(a) is then described as 



LX : \xo, X1.X2. . . . ,x. 



\xo,Xi,X2, ■ ■ ■ ,Xn) (Bl) 

where the phase (pj+i in the exponent depends on the in- 
teraction time and the interaction strength between two 
atoms at the lattice site j + 1, and the addition is per- 
formed modulo 2. Note that two neighboring atoms at 
the sites j and -f 1 contribute to the exponent if and 



only if Xj = and 



--j+i 



= 1. The variables Xj can only 



take on the values and 1. In all examples we discuss 
here, ipj ^ ip = constant and is the same for all lattice 
sites. 



The operation (Bl) defines a generalized phase gate 
that acts on a group of n neighboring atoms via shifting 
the lattice across one lattice site. It can easily be seen 
that, for example when combined with 7r/2-pulses as in 
Fig. 10, LX produces the entangled states ( ^3| ) and ( |3^ ) 
for n — 2 and n — 3. 

In two dimensional lattices, the logical variables Xki are 
labeled by two indices, where k is the horizontal index 
and I the vertical index. The phase gates corresponding 
to horizontal and vertical lattice shifts are then defined 
as 



mod2)a;j + i,,ipj + i_i 



LX\{xki})^e 
LY\{xki}) ^ e-'^^'"'^+''''°''^^'''-+''^'-+'\{xki}) (B2) 



as an obvious generalization of (Bl). 

It is clear that the operations can be further gener- 
alized to lattice shifts across an arbitrary number of 
lattice sites and along arbitrary directions. There are 
interesting topological questions in this general situa- 
tion. For the present discussion, however, the gates 
LX, LY as defined in (B2) are sufhcient and we will 
set ipk^i = TT. Apart from these gates, we will only 
need single-particle operations, in particular the Pauli- 
operators <Jxj,cryj,azj and the Hadamard transforma- 



)/V2 applied to an 



tion (7r/2 pulse) Hj 
atom with index j. 

Consider now a configuration of 3 x 3 atoms as in 
Fig. ^ where the central atom is in the unknown state 
\ip) = q;|0)-|-/?|1) while all surrounding atoms are initially 
in the state |0). Let us first look at the special case when 
lip) = |0), that is, the central atom is in the state |0) as 
well. If we apply a 7r/2 pulse to each atom of the block 
and then the operation LX, we obtain a tensor product 
of three GHZ states where each row of the block is in the 
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same state (0 + 1)0 (0-1)- (0-1)1 (0+1). (For notational 
brevity, we suppress the bracket notation in the follow- 
ing and identify = |0) and 1 = |1)). This state can 
be transformed to the form 000— 111 by applying Hi to 
the first atom and H^az^s to the third atom of each row. 
The operation LX, supplemented by one-qubit rotations, 
produces thus one of the code words in (^. 

To realize a quantum memory, an unknown state '0 = 
aO + /31 of the central atom is to be encoded into an 
entangled 9-bit state as in (^7|). Let us write the initial 
(unencoded) state of the block in the form 



for each error the corresponding syndrome and the state 
of the central qubit: 



Ibare) = O1O2O3 OiiaOs + PWOe OrOgOg 



(B3) 



where the first, second, and third triplet refers to the 
upper, center, and lower row of the block in Fig. |l^. 
To encode into a corresponding superposition of both 
codewords, lattice movements in both horizontal and in 
vertical direction are required. In detail, the encoding 
operation is given by 

ENC = Hie o LX o H^^q o LY o axisedHiaievg o LX o H"" 

(B4) 

whose essential part is a sequence of three lattice move- 
ments, horizontal - vertical - horizontal, with 
certain 1-bit unitary transformations in between. In the 
notation used here, denotes a Hadamard transform 
applied to each of the 8 syndrome atoms, whereas Hijk... 
and az-ijk... are single-qubit rotations applied to th e se - 
lected atoms i, j, fc, . . ., only. Applied to the state (|B3| ), 
ENC produces 

ENC Ibare) = a(000 - 111)(001 - 110)(000 + 111) 
/3(000 + 111)(100 + 011)(000 - 111) 
= «0l+/31l. (B5) 

The codewords Ol and 1l are equivalent to the Shor code 
(p6|), as we shall see presently. 



The decoding operation is given by the inverse of (B4), 

DEC = iJ^ o LX o ffi346790-2;:369 ° LY o iJ45g o LX o H^q 

(B6) 

involving the same lattice movements, but the 1-bit oper- 
ations carried out in reverse order. To see explicitly how 
one can correct an error occurring on one of the qubits 
j = 1, 2, . . . 9 , we apply the error operators axj, (^y.j, or 
CTzj- to the encoded state aOL + /31l- Then we apply the 
decoding operation DEC and measure the state of the 
syndrome atoms. The code Ol and 1l is error correcting 
if every possible error is mapped into a syndrome state 
different from O1O2O3 O4O6 OyOgOg, and for each syndrome 
we can tell which unitary transformation has to be ap- 
plied to the central qubit to restore it to its original state 
[it is not necessary that all errors are mapped to mutu- 
ally orthogonal subspaces [0]. The following table gives 



error 


syndrome 


central qubit 


none 


Ann nn nnn 
UOOUU UUU 


aO 


+ P 1 




no 00 000 


aO 


-13 1 


(yx.2 


1 n 1 A A AAA 

lUi UU UUU 


aO 


-/3 1 


fx, 3 


A1 1 AA AAA 

Oil 00 000 


aO 


-P 1 


fa:, 4 


01010010 


aO 


+ (3 1 


fa;, 5 


UUU ii UUU 


aO 


-P 1 




UiU Ui UiU 


aO 


+ (3 1 


l^xj 


000 00110 


aO 


-13 1 




Ann nn i m 
000 00 101 


aO 


-13 1 




UUU UUUii 


aO 


-13 1 




100 00000 


aO 


-13 1 




111 AA AAA 

111 00 UUU 


aO 


-13 1 




A A 1 A A AAA 
UUl UU UUU 


aO 


-f3 1 


<Jy,4 


000 01000 


al 


- P 


0-y,5 


A1 A AA A1 A 
UiU UUUlU 


al 


-(3 


<Jy,6 


AAA 1 A AAA 

000 10000 


al 


- (3 


^V,7 


000 00100 


aO 


-(3 1 




AAA AA 111 

000 00 111 


aO 


-(3 1 


(JyS 


AAA AA AA1 
000 00 UUl 


aO 


-(3 1 


(^z.l 


010 00000 


aO 


+ (3 1 


0'z,2 


A1 A AA AAA 

010 00 UUU 


aO 


+ (3 1 




A1 A AA AAA 

UlUUU UUU 


rv/O 


J- R 1 




01011010 


al 


-/3 




01011010 


al 


+ /3 


(^z,6 


01011010 


al 


-(30 


(^z,7 


00000010 


aO 


+ (3 1 


^zfi 


00000010 


aO 


+ (3 1 


(^z,9 


00000010 


aO 


+ (3 1 



In Fig. |I3|, the syndrome atoms visually encircle the 
unknown qubit that is to be protected. If any of the 9 
qubit suffers a spin flip, a phase flip, or both, the error 
can be detected by measuring the state of the syndrome 
atoms after the decoding operation has been applied to 
the group. This could be done by a fluorescence mea- 
surement where atoms in state 1 and correspond to 
"bright" and "dark" , respectively. For example, accord- 
ing to above table, the pattern 



1 V' 1 


tells us that a spin flip has occurred in the central atom, 
whereas 
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V' 

1 1 

reveals a spin flip in the left atom of the lower row, and 

1 

1 V' 1 
1 

corresponds to a phase flip in any of the atoms of the 
central row. In any case, the state i}}' of the central qubit 
after the detection of an error is related to the initial state 
■0 via a (known) unitary operation U: ip' = U^'^ which 
can be obtained from the third column of the syndrome 
table given above. 

The fact that the encoding operation involves only 3 
lattice movements provides a specific example of a "par- 
allelization of a quantum circuit" We have not 

proven that 3 is really the minimum number of entan- 
glement operations needed; there might be still faster 
sequences. The original Shor code can be recovered from 
this code by applying an additional vertical lattice shift, 
LY, and certain 1-bit rotations. 
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